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Abstract 

We present a non-perturbative renormalization of the bound state problem of n bosons 
interacting with finitely many Dirac delta interactions on two and three dimensional 
Riemannian manifolds using the heat kernel. We formulate the problem in terms of a 
new operator called the principal or characteristic operator <&{E). In order to investigate 
the problem in more detail, we then restrict the problem to one particle sector. The 
lower bound of the ground state energy is found for general class of manifolds, e.g., for 
compact and Cartan-Hadamard manifolds. The estimate of the bound state energies in the 
tunneling regime is calculated by perturbation theory. Non-degeneracy and uniqueness of 
the ground state is proven by Perron-Frobenius theorem. Moreover, the pointwise bounds 
on the wave function is given and all these results are consistent with the one given in 
standard quantum mechanics. Renormalization procedure does not lead to any radical 
change in these cases. Finally, renormalization group equations are derived and the j3 
function is exactly calculated. This work is a natural continuation of our previous work 
based on a novel approach to the renormalization of point interactions, developed by S. 
G. Rajeev. 

1 Introduction 

The studies of Dirac-delta interactions in quantum mechanics (which are also called zero 
range, contact or point interactions, or Fermi pseudopotentials in the literature) date 
back to the work of Kronig and Penney pQ who introduced the periodic delta interac- 
tions describing the non-relativistic electrons moving in a one dimensional fixed crystal 
lattice. The historical development of point interactions has been given extensively in the 
monograph [2]. Although it was Thomas [3] who pointed out that the problem of point 
interactions in three dimensions could not be physically acceptable due to the ultraviolet 
divergences, Thorn [4] realized that we did not have to abandon these interactions and 
physical results could be obtained after regularization and renormalization, well known 
procedures in quantum field theory. 
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The motivation of studying point interactions in two and three dimensions is based on 
trying to understand the concept of renormalization in a simple context rather than field 
theory. There are large amount of works on the renormalization of point interactions in the 
literature from several point of views [3] - [11]: Regularization schemes can be performed 
either in coordinate space [8] or in momentum space [H [TOj, [TH [121 CC2 O IS] • 
A single point interaction in two dimensional flat space is also an instructive example 
of dimensional transmutation [H [9[ fT6[ I17j . that is, the original Hamiltonian does not 
contain any intrinsic energy scale due to the dimensionless coupling constant in natural 
units, but a new parameter specifying the bound state energy, is introduced after the 
renormalization procedure, which then fixes the energy scale of the system and this is 
called dimensional transmutation. This implies a violation of SO(2, 1) symmetry of the 
scale invariant potential, so it is one of the simplest examples of anomaly or quantum 
mechanical symmetry breaking [10] . Furthermore, renormalization group equations of 
point interactions have been discussed in [5lll5[IT8] and the (3 function has been calculated 
exactly so that the theory has been found as asymptotically free in two dimensions. 

Mathematically rigorous treatments of point interactions are given in the context of 
self-adjoint extension theory [19] and a detailed exposition of this subject has been dis- 
cussed in the monograph [2]. Dirac delta interaction is considered as a self-adjoint ex- 
tension of a formally Hermitian free Hamiltonian on a space with a removed point. The 
result of this method is identical to that of the renormalization method if a certain relation 
between the parameter in the extension and the renormalized (or bare) coupling constant 
is satisfied [10] . Self-adjoint extension method can be also analyzed within the Green's 
function method [2T71 f2T] . 

Many body version of this problem on R 2 and R 3 is known as the formal non-relativistic 
limit of the Ac/) 4 scalar field theory in (2+1) and (3+1) dimensions. All these are extensively 
discussed first in the unpublished thesis of J. Hoppe [22] and later from a new perspective 
in [231 124], S. G. Rajeev [23] introduced a new non-perturbative renormalization method 
for point interactions which can be applied to several many body theories: quantum 
mechanics with point interactions, fermionic and bosonic quantum fields interacting with 
a point source, many body problems with point interactions, and non-relativistic field 
theory with polynomial interactions. One of the main advantages of this approach is that 
all the information about the spectrum of the model is described by an explicit formula 
instead of imposing the boundary conditions on the operator as in the case of self-adjoint 
extension theory. This method is also particularly useful for dealing with the bound state 
problems because of its non-perturbative nature. We are not going to review the ideas 
developed in there. Instead, we suggest the reader to read through the paper [23] to make 
the reading of this paper easier. 

Following the original ideas developed in [23J, we previously considered the bound state 
problem for iV point interactions in two and three dimensional Riemannian manifolds [25]. 
The lower bounds on the ground state energy are found for few special manifolds §> 2 , H 2 
and H 3 [25]. The construction of the relativistic extension of this model in two dimen- 
sions is also possible [26]. We also applied the same method based on the heat kernel to 
non-relativistic Lee model [27] and its relativistic version has been constructed later on 
[28j . Our primary motivation comes from the question, how the renormalization method 
for the point interactions in quantum mechanics should be performed non-perturbatively 
on Riemannian manifolds, hoping that we can extend our understanding to the realm of 
quantum field theory. It is worth pointing out that this problem on two dimensional Rie- 
mannian manifolds also displays a kind of dimensional transmutation, where new energy 
scales different from the intrinsic energy scales of the system appear after the renormal- 
ization 1251 . 



We organize this paper as follows. In section [21 we construct the bound state problem 
for n bosons living in two and three dimensional Riemannian manifolds interacting with 
N external Dirac delta interactions. This construction is motivated by the work [23] in 
which the many body version of this problem is renormalized non-perturbatively in flat 
spaces. We then restrict the problem to n = 1, and find the finite formulation of the 
problem in terms of a new operator, which is called principal operator (or characteristic 
operator) [23| and the result is consistent with the one that we have obtained in }25j . 
In section [3l we find the wave function for the bound states and show that ground state 
energy of iV + 1 center case is smaller than the N center case using Cauchy interlacing 
theorem. Then, we calculate the bound state energies in the tunnelling regime using a 
version of perturbation theory. The results for our own purposes on the upper and the 
lower bounds of heat kernel in the mathematics literature is given shortly for compact and 
Cartan-Hadamard manifolds in subsection 14. II and in 14.2} respectively. Section [5] presents 
pointwise bounds on the wave function using the upper bounds of the heat kernel given in 
the previous section and it is shown that the result is consistent with the classical result 
found in the standard quantum mechanics. Section [6] establishes the lower bound of the 
ground state energy for more general class of Riemannian manifolds. The proof of the 
lower bound for the ground state energy has the same spirit with the one given in our 
previous work |25] but the proof given here is generalized to a rich class of manifolds, 
such as compact and Cartan-Hadamard manifolds. Then, in section [71 non-degeneracy 
and positivity of the ground state is proven with the help of Perron-Frobenius theorem. 
Finally, we proceed with the study of the renormalization group equations and the j3 
function is calculated exactly. Some important properties and asymptotic expansion of 
the heat kernel is summarized in the appendix A and the proof of the existence of the 
Hamiltonian in two dimensions is explicitly given in the appendix B. 



2 Construction of the Model 

We consider n nonrelativistic bosons living in a two or three dimensional (D = 2, 3) 
Riemannian manifold and they interact with N external attractive Dirac delta potentials. 
In the second quantized language, Hamiltonian of the system is 

JM 

where d®x = \/det g d D x is the D dimensional volume element, is Laplace-Beltrami 
operator in the local coordinates x = (x 1 , x 2 , . . . , x D ) 




and <jy g (x), <p g (x) is defined as the bosonic creation-annihilation operators on the Rieman- 
nian manifold with metric structure g. Here a% denotes the location of i-th external dirac 
delta potential on the manifold and \ € M + is the strength of the delta interaction at Oj, 
so called coupling constant. 

It is easy to show that the number of bosons f M d^x <Pg(x) 4> g {x) is conserved. The 
lowest eigenvalue of the Hamiltonian H of our problem in a sector with fixed number of 
bosons is either zero or negative infinite, so the energy of the ground state is not bounded 



4>\{x) 



2m 



N 

^2Xi4>l(x)6g(x,ai)(/)g{x) 
1=1 



(1) 



3 



from below |25j : E — > — oo. The first step we must do is to regularize the model. The 
natural regularization of Hamiltonian can be chosen as 



N 

H e = H Q - Ai(e) / dfxdfy K e/2 (x,ai]g)K £ / 2 (y,ai;g)<f>l(y)<j)g(x) , (3) 
i=i Jm2 

where Hq is the free Hamiltonian. K e (x, y; g) is the heat kernel defined on the Riemannian 
manifold (M,g) and it converges to Dirac delta function 5 g (x,y) as e — > + . In this limit, 
one can see that we recover the original Hamiltonian we are interested in. We write the 
heat kernel as K e (x, y; g) throughout the paper in order to specify which metric structure 
it is associated with. Some essential properties of the heat kernel on Riemannian manifolds 
that we have used in this paper are given in the appendix A. 

Now, we will consider the resolvent of the regularized Hamiltonian in a Fock space 
formalism with arbitrary number of bosons. Following the same methodology developed 
for the model in the plane [23], we shall extend the bosonic Fock space B that we have 
started with, to B = B © B ® by defining new creation and annihilation operators 
at the locations of the dirac delta interactions. These are called angels, which is first 
introduced in [23]. The angel states allow us rewrite the model in such a way that the 
coupling constant appears additively rather than multiplicatively. As a result, we can 
renormalize the model nonperturbatively by simply normal ordering. We assume that 
the angel operators obey the orthofermionic algebra [29\ defined by the following product 
relations (not with commutators): 

N 



:« x] + &ij x l Xk = i<% , Xi xj = o = x 1 x] . ( 4 ) 

k=l 



where 1 is the identity operator and k = 1,2,... ,N. It is more convenient for our 
purposes to write the angel algebra in terms of projection operators: 

Xi x) = SijUo, Xi Xj = = x\ x) , (5) 

where 

N 

n i = ^xix/t, no^i-nx, (6) 

k=l 

are the projection operators onto the 1-angel and no-angel states, respectively. Now we 
define the augmented regularized Hamiltonian H e on B 



N , N 1 

H e = H U ^ V/ dgX K e / 2 (x,ai\ g)(f> g {x)x\ + h.c + J2 i-yrxl Xi ■ (7) 
i=i J M J i=1 A ^e) 



If we split the Hilbert space according to the angel number, the corresponding operator 
H € — EHq can be written in the following matrix form: 

with a : B -> B, bt : B © -> B, d £ : B <g> -> B ® C N . Here, 

N 



Ho-E, de = J2^ 7) xhi 

i=i * w 
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N 

b\ = y2 / dg x K e / 2 {x,ai;g)4>\{x)xi ■ 



(9) 



Then, one can construct the augmented regularized resolvent 

1 



Re{E) 



H e - Eli 







a e 01 



(10) 



One can find a e ,/3 e ,5 e in terms of a,b e ,d e by direct computation. This could be done 
apparently different but equivalent ways and the formulas were obtained in the appendix 
of [23] 

Re(E) . (11) 



a-b\ d' 1 b e 



H f — E 



This means that R t {E) projected to B is just the resolvent of the operator H t . We have 
also another formula for a e 



a 1 + a 1 bt d e -b e a 1 b\ 



b f a 



or 



where 



Re(E) 



1 1 

+ 



Hn — E Hn — E 



b\^ e {E)-\ 



Hn — E 



(12) 
(13) 



N 



i=l 

Ay) 



Hn — E 



~ Yl / dfxdfy K e/2 (x,a i ;g)K e/2 (y,a j ;g) 

4>\{x)x\xj ■ 



(14) 



The operator & t {E) is called the regularized principal operator (or regularized character- 
istic operator [23]). Note that writing the resolvent of H e in this way allows us to write the 
coupling constant additively. The renormalization procedure can then be done if we can 
separate the singular part of the operator <J> e . We will see that this is possible by normal 
ordering of the operators in the principal operator. By using eigenfunction expansions 
(|175p and (|177p of the operators 4> g {x), (jy g (x) and that of the heat kernel (or their analogs 
for non-compact manifolds), one can shift the operator (j)\(x) in (|14p to the left 



Ha-E 



•to 



J dfx' ^I(X') jH l e"* <*«>-*> Kt(X: X >. g) 



(15) 



and one can also shift the operator 4> g (x) to the right using a similar equation. Then, the 
normal ordered principal operator can be written by using the properties of heat kernel 
and separating the i = j term from the sum 

N 1 N f f°° dt 

~[ M e ) ~Z 1 Jm 2 Jo h 



i=l Jo 



£/° ^ (e+t) (a i ,a i ; 5 )e-l^) x t Xj 



(16) 



i=l 

i+3 



Due to the singular behavior of the diagonal heat kernel (|18ip near t = 0, we expect that 
the third term is divergent as e — > + . Therefore, if we choose the coupling constant 



dt -K t (a i ,a i ;g)e-T^ , (17) 



M e ) Je ft 

where — //? corresponds to experimentally measured bound state energy of the individual 
ith. Dirac delta center, we find the principal operator after taking the limit e — > + 

N 

N 



$(£) = lim $ e (E) = V / -K t K, a 4 ; 5) 
_ Z) / IT/ df«dfy^t(ai,a:;i/)iirt(o i ,!/;p)4(x)e-K^ -^^(y)xIxj 

-E/ l^^;^-^ -^^- (is) 
i=i ^° 



This can be written in a more compact way $>(E) = X]fj=i ^ij(E) x\ Xj> where $>ij(E) 
can be read from (|18|) . Once we have a proper definition of the principal operator, the 
divergence is completely removed since the spectrum of the problem can be found from 
the resolvent. We are now in a position to get the full resolvent of our problem in terms 
of the principal operator 

1 N N 

R(E) = to R,(E) = + — J>jM Xt •"'(*) 

fc = l i = l 

11^ 1 

= + ^ S •5 1 W* W 7^ • (•») 

i,3=l 

where we have used ^~ 1 (E) = X2fj=i ^H-^) xl Xj an d the algebra of angel operators © 
with the fact that R(E) : B — > B. We note that the principal operator can be extended 
to its largest domain of definition in the complex energy plane by analytic continuation. 

Since R(E) : B — > B, we can consider the resolvent kernel between n bosons with no 
angel states. Up to here, we have generalized the construction of the problem given in [25] 
to the many boson cases. From now on, we shall study one boson problem for simplicity, 
that is, 

/ dfxi(>(x)4>l(x)\o)®\n), (20) 

J M 

where |0) is the vacuum for the angel state and ip{x) is the wave function for the boson. 
Then, the resolvent kernel corresponding to this state satisfies the following equation after 
a straightforward calculation 

N 

R(x,y\E) = R (x,y\E) + £ R (x,a i \E)^>r j 1 (E)R (a j ,y\E) , (21) 

where Ro(x, y\E) is the free resolvent kernel which can be written in terms of heat kernel 
(I180p . Here &ij(E) can be analytically continued to its largest set in the entire complex 
plane so E should be considered as a complex variable. The equation (|2ip gives the 
relation between the resolvent defined on an infinite dimensional space and the principal 



matrix defined on a finite dimensional space. Such formulae were extensively discussed in 
problems associated with self-adjoint extensions of operators, notably by Krein and his 
school, and also for such singular interactions in flat spaces [21 [TO] 121] . Therefore, our 
problem can be considered also a kind of self-adjoint extension of the free Hamiltonian. We 
will come back to this point at the end of the section [5j The resolvent essentially includes 
all the information about the spectrum. We will restrict ourselves only to the bound state 
problem since the scattering problem requires a deeper analysis. The discrete spectrum 
of the Hamiltonian is the set of numbers E at which the resolvent does not exist and 
continuous spectrum corresponds to the unbounded resolvent. The poles corresponding 
to bound states must be due to the <1> -1 (.E7) which can be seen from (|2ip . In other words, 
the roots of 

=0, (22) 

determine the bound state spectrum of the model. Since $>(E) : B ® C N —IB® C N , let 
us try to consider |^) as a direct product of no-boson with one angel state: 

N 

l*> = \0)®^A k \e k ) , (23) 
k=l 

where \ek) = xt|^) is a set of complete orthonormal basis for C^. Then, the equation 
(1221) yields 



V / —K t ( ai ,ai;g) (e-i^ - ei E ) Ai\ ei ) - J2 / -K t (a h a j ;g)eTi E A j \e i )=Q,(2&) 
l=1 Jo n v > fr!J n 

and the result can be written as a matrix equation 

N 

^2^ ij (E)A j = Q, (25) 



where 



J ^-K t (ai,ai]g)(e "rrf -eh E ^j if i = j 

° f°° dt t 
-J -^K t (ai,aj;g)e^ E tfi^j- 



(26) 



Thanks to the symmetry property of the heat kernel (|173h , the principal matrix is Hermi- 
tian for real values of energy and the explicit form of it in terms of heat kernel has been 
first obtained by a different method in our previous work [25]. Many important aspects 
of the model can be understood by working out the principal matrix as we will see in the 
following sections. 

It is well known that the same problem with a single delta potential in flat spaces is a 
good example of a dimensional transmutation in quantum mechanics. Our problem in D = 
2 realizes a generalized dimensional transmutation [17] : In our case, the coupling constants 
Xi have the same dimension as ^- by dimensional analysis. In contrast to the flat case, we 
also have intrinsic scales coming from the geometry of the space, such as curvature and 
the geodesic distance between centers dij. However, after the renormalization procedure, 
we obtain a set of new dimensional parameters /x? from the relation ()1T[) . Hence, the 
energy is not determined by naive dimensional analysis. However, in the case of single 
delta attractor for the flat case there is no combination of dimensional parameters to 
come up with an energy scale, whereas in the case of a manifold we have geometric length 
scales which already may define an energy scale. The dimensional transmutation is most 
striking in cases where there is no intrinsic energy scale. 
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3 Interlacing Theorem and Perturbation Theory 



A mathematically satisfactory calculation of the wave function should proceed from the 
resolvent equation. Since the eigenvalues are isolated we can find the projection operator 
to the subspace corresponding to this eigenvalue by a contour integral |30j : 



(x\F k \y) =M*W k (y) 



t^t jf dzR(x,y\z), 



(27) 



where r& is a small contour enclosing the isolated eigenvalue —vi. We note that the free 
Green's functions Ro(x,y\z) will not contain any poles on the negative real axis, so all 
the poles on the negative real axis will come from the poles of inverse principal matrix 
^~ 1 (z). For simplicity we can assume that the eigenvalues are nondegenerate and let us 
denote the k th eigenvalue of the principal matrix as uj k so eigenvalue problem is 



N 



Since the principal matrix is Hermitian on the real line and 

*J,(z) = <l>y(z*), 



(28) 



(29) 



on the complex plane, there exists a holomorphic family of projection operators on the 
complex plane [31], so that we can apply the spectral theorem for the principal matrix 

$y(z) = J> fe (z)P fc (% , (30) 
k 

here Pf.( z )ij = A k *(z)A k (z), Af(z) is the normalized eigenvector corresponding to the 
eigenvalue u k (z). Similarly, we can write the spectral resolution of the inverse principal 
matrix, 



The residue can then be found 

lim Ro(x,ai\z)(z + v%)$rMz)Ro(a,j, y\z) 



R (x,ai\ - v\) 



' du) k (z) 




dz 





\(-Vk)ij R o{aj,y\ - vl) ■ 



(32) 



Now we will look at the variations of the eigenvalues of as we change the parameters v. 
Using 

(33) 



= [A k (-y 2 )^{-v 2 )A k (-v 2 )) 



and as a consequence of Feynman-Hellman theorem [32] in the non-degenerate case, we 
have 



du) k (—v' 2 
dv 



A k (-u 



2 , <9$(-^ 2 ; 



N 



2 .d^(-u 2 ) 



A k (-u 2 ) 



Oi 



A k (-u 2 ) 



(34) 



Taking the derivative of the principal matrix with respect to v from ((26 



Oi 



v=v k 



dt 



( 2v h t 



K t (ai,aj;g)e ft 



to h \ h 

and inserting the equation ([35]) into the equation (fM|h we obtain 
doj k (z) _ 1 du k (-u 2 ) 



(35) 



1 

2^ 



dz 

N 



2v di 



E A i*(-» 2 ) 



dt /2z/i\ 



o h \ h J 



)K t (a i ,a j ;g)e * A k A-v 2 ) . 



(36) 



If we evaluate (|36|) at z = — or at i/ = Uk, it yields 



^tK.Oji^e s-Aj(-i^). (37) 



Note that the integral is finite in two and three dimensions due to upper bounds on the 
heat kernel. If we combine all these results, we get 



^k{x)ilf k {y) 



1 

2m 



-^(27ri)-Ro(a:,ai| - 



xK t (ai,aj;g)e ft Aj(-v k ) 



(38) 



: t (ai,aj\g)e * Aj(-v%) 



Then, we can directly read off the bound state wave function from the equation above: 

xy e—fT22Ai(-ul)K t (ai,x;g) , 



(39) 



i=l 



from which one can easily see that ipk( x ) is finite except at x = at. This is exactly the 
same result obtained by a different method in [25j . Note that we have written A k (—u 2 ) 
as Ai{—v 2 ) for simplicity. Incidentally the equation (f37|) implies an interesting result for 
the variation of eigenvalues, 



duj k {—v 



dv 

N 



E A *^[% (it) 

^2 A *i(~^k) -j- \~YJ J d 9 x K t/2{a i ,x;g)K t/2 (x,a j ;g)e~~^A j (-^l) 

/" ^ 2 

/ d g x y2\ K t/2(ai,x;g)Ai(-ul)\ , (4 



dt (2v k t\ _*£ 
o h \ H 



where we have used the fact that v k G ffi + for all A; and properties of heat kernel (|172[) 
and (|173p and the order of the integration and the finite sum can be interchanged. We 
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can easily see that the above equation (|40p is strictly positive due to the positivity of heat 
kernel, so 

du k {-v 2 ) 



Oi 



> . (41) 



Energy eigenvalues E = —v are obtained from the zeros of the eigenvalues of the principal 
matrix, that is, uj k (—v k ) = 0, and there is a unique solution for each uj k {—v 2 ). We also 
know that for sufficiently small values of u, the matrix $>(— v 2 ) becomes negative, hence 
no zeros exists beyond some critical point. 

With this result in mind, let us see what can be said about the comparison of the 
energy eigenvalues for the different number of delta centers. In order to see this, we need 
the Cauchy interlacing theorem in mathematics literature [33 j . which states that if we 
delete the last row and column of an Hermitian (N + 1) x (N + 1) matrix, the eigenvalues 
of the original matrix is interlaced by the eigenvalues of the new matrix, i.e. if oj {— v 2 ) < 
lo 2 (—v 2 ) < ■ ■ ■ < lo n+1 (— v 2 ) lists the eigenvalues of the original (N + 1) x (N + 1) matrix 
and if u 2 ) < Cj 2 (—i> 2 ) < ■ ■ ■ < oj N (—u 2 ) lists the eigenvalues of the reduced matrix 
(any N x N principal submatrix of the (N + 1) x (N + 1) matrix), then we have 

oj\-v 2 ) < Cj\-v 2 ) < oj 2 {-u 2 ) < Cj 2 {-u 2 ) <■■■< Cb N (-u 2 ) < u N+l (-u 2 ) . (42) 

We assume that (iV + 1) x (iV + 1) matrix in the above-mentioned theorem is the principal 
matrix <fr N+1 (—is 2 ) corresponding to a certain arrangement of N + 1 delta potentials. 
If we now delete the last row and the column it means that we remove the (N + l)th 
delta center from the system. The bound state problem of N centers corresponds to zero 
eigenvalue of the principal matrix $ , and let us denote that bound state energy as E: 

u k (-9* k 2 ) = 0, E k = -Df, (43) 

if it exists. By the Cauchy interlacing result, we then expect the following inequality 

• • • < u k (-9* k 2 ) < Cj k (-9* k 2 ) = < o/ +1 (-^ 2 ) < • • • , (44) 



From the positivity of the derivative of the eigenvalues with respect to the argument (|41 j) . 
the eigenvalues oj are monotonically increasing functions. Hence, in order to get a zero 
root of Lo k (—i> k 2 ), we should increase i> k to a higher value v k . As a result, Lo k (—u k 2 ) = 
if u k 2 > u k 2 , i.e. E k = —v k 2 < E k = — u^ 2 . Thus, the energies also interlace in the same 
manner-this is a nonlinear analog of Sturm's comparison theorem of the eigenvalues. 
Moreover, E^ r +1 < E^ r < E x gr = — [i k , that is, the ground state is always negative and 
approaches the bound state energy for the one delta center as iV gets smaller. We can 
also generalize these results to the degenerate cases but the proof is more cumbersome. 

It is worth noting that we do not have to solve the energy eigenvalues for the bound 
state while performing our non-perturbative renormalization method. In other words, 
although our renormalization method makes the problem well defined and finite, the 
energy eigenvalues must be found after this finite formulation has been constructed. If we 
cannot solve the problem exactly after the renormalization, we must apply the standard 
approximation methods, such as perturbation theory and variational techniques. An 
interesting estimate for our problem can be given by perturbation theory. For simplicity, 
we assume that all binding energies — [i 2 s are different and the magnitude of the minimum 
binding energy of the k-th singular potential is much larger than the correlation energy 
between the k-th fixed center and the Z-th center, that is to say, we assume 2m <$(a fc a t ) 
[i k min (on a noncompact manifold we may assume that the geodesic distance between the 
centers is large) . This assumption makes the off-diagonal elements of the principal matrix 
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much smaller than its diagonal elements. For this reason, let us separate the principal 
matrix for E k = —v\ as the sum of a diagonal matrix and an off-diagonal matrix, which 
is very small compared to the diagonal part: 



(45) 



Since <&(— is Hermitian, we can apply standard perturbation techniques to our problem. 



The eigenvalue problem for the principal matrix we wish to solve is given in (J28J). We 
again suppose that there is no degeneracy for simplicity. The energy eigenvalue changes 
to E k = E { k 0) + 5E k or 



v k = vf + Sv k ■ 



(46) 



From the fundamental idea of the perturbation theory in finite dimensional spaces, we 
have the following expansions for the eigenvalues and eigenvectors: 



UJ 



A fc(0) + A fc(1) + A fc(2) + -- - , 



(47) 



and the solution to the related unperturbed eigenvalue problem 

N 



£[M-^4 (0) (-^ 2 ) 



fc(0) 



(48) 



is given by 



CO 



fc(0) 



(-4) 







jK t (a k ,a k] g)[e- t ^ n 



-tv 



!/»] . 



(49) 



Then, the energy eigenvalues can easily be found from the condition uo k ^\— ui) = 0: 



E 



(o) 



-t4 



or 



,(°) 



(50) 



and eigenvectors are 



A k(0) 



(51) 



where 1 is located in the A;th position of the column and other elements of it are zero or 



we can write A 



fc(0) 



S^. Here el s form a complete orthonormal set of basis. 



N 

E 

i=i 



(52) 



We must emphasize that there is no distinction between upper and lower indices for our 
purposes. Once we have found the solution of the diagonal part of the principal matrix 
or unperturbed eigenvalue problem, we can perturbatively solve the whole problem. The 
standard perturbation theory gives us the first and second order eigenvalues: 



N 



■UJ 



fc(l) 



1 1 



N 



id 



fc(2) 



= E 



N 



- w*«J)(-i^)-w'(o)(-^) 



E 



*»(-^)*H(-ffc) 



(53) 



respectively. Hence the energy eigenvalues of the whole problem can be determined from 
c/(-/4 + 5E k ) = u> k ®(-n 2 k + 6E k )+cj k( - 2 \-f4 + 5E k ) + . . . = , (54) 
and uj k( - \— z/|) and <3?fcz(— ff) for ^ I can be expanded around v k = [i k 



<S> kl (-f4 + 5E k ) 



dv k 



5v k + 0{5 2 u k ) 



Ji* + 0{8 2 u k ) , (55) 



where we have used the fact uo k( -°\— fi 2 .) = 0. If we substitute ([55]) into (f54"|) and ([53]) . and 
use Feynman-Hellman theorem (|34|) . we obtain 



JV 



l^k 



$ k i(-f4)$i k (-i4) 



i + 



dv k 



Vk=P-k 



dv k 



du k 



Vk=Vk 
1 



Vk=Vk 



8u k 

+ 0(5 2 v k ) . (56) 



If we also expand the last factor in the powers of 5u k and ignore the second order terms 
and combine the terms using the symmetry property of principal matrix, we find 



dv k 



N 

^fc=£»fc 1=1 
l^k 



dv k 



Vk=Vk 



N 

£ 

i=i 

l^k 



dv k 



®ki{-V 2 k )<5>ik{-H 2 k ) 



N 



3> k l(-vl) d^ik{-y 2 k ) 



Vk=Vk 



Vk=Vk 



Su k 



+ 0{5 2 v k ) . 



(57) 



Ignoring the second and third terms on the left hand side of the equality due to the fact 
that § kk (— v k ) 3> \& k i(— v k )\, we get the change in v k 



bv k 



d<5>kk{-v 2 ) 



dv k 



\- l N 

E 

Vk=P-kJ 1=1 
l^k 



+ 0(5 2 v k ), 



(58) 
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so the change in the energy is dE^ ~ —2^k8u^ + 0{5 v^). Let us now consider how the 
bound state energy changes in the tunneling regime for our problem in which 2~jr "C 

fJ'kmin- We mus t nrs t calculate asymptotic behavior of the off-diagonal element of the 
principal matrix in this regime. In order to see this, we make the scaling transformation 
t = u/B, where B = h/2md 2 j and use the scaling property of the heat kernel (I179j) . In 
two dimensions, we have 

®ij{~l4) = ~ j Q J- R u(ai, af, Bg)e~^B . (59) 

In tunnelling regime, the most significant contribution to the integral comes from the 
region u = due to the fact that integrand is suppressed by the exponential term for 
large values of u. Hence we can use the short time asymptotic of the heat kernel given in 
(1 1 841) . The result is an integral representation of the modified Bessel function of the third 
kind 



where we have used dij — > B l / 2 dij and ^ — > ^1/B 1 ^ (for two dimensions) under the 
scaling transformation g — > Bg. The asymptotic expansion of Kq(x) for large values of x 

Ko(x) ~ J^e-' (61) 



leads to 



Here, large values of x corresponds to tunneling regime in our problem. For three dimen- 
sional case, the idea is the same and the result would be 



-yl^S*^')*-^. (63) 
[l mh {Aixh/2mf/ 2 y ' 

Therefore, in the tunnelling regime, we can find the change in the bound state energy 
8Ek in the presence of other delta interactions by substituting (|62p and (|63p into (|58p . In 
agreement with the naive expectation in the standard quantum mechanics, we show that 
the bound state energy in the tunnelling regime gets exponentially smaller with increasing 
distance between the centers. 

Before we prove the lower bound for the ground state energy and find pointwise bounds 
on the wave functions, the results in mathematics literature for the upper and lower bounds 
of the heat kernel for compact and Cartan-Hadamard manifolds will be briefly given in 
the next section. 



4 Upper and Lower Bound Estimates of the Heat 
Kernel 

4.1 Compact Manifolds: 

1) Upper Bound of Diagonal Heat Kernel: The following result is a simplified version of 
the corollary 3.6 given in [35], which assumes that some geometrical conditions must hold 

13 



on the boundary. The global upper bound estimate of the diagonal heat kernel in [35] 
includes whole boundary information via an explicitly calculable strictly positive constant 
A = A(diam(.M), H, R, K, V (Ai)), where diam(.M) is the diameter of the manifold, and 
K is the lower bound on the Ricci curvature, and also H and R are parameters related 
to boundary conditions. We then state the following corollary by safely removing the 
boundary effects in corollary 3.6 in |35| since A is strictly positive: 

Let Ai be a compact manifold. Suppose that the Ricci curvature of Ai satisfies 
Ric M > -K, K > 0. Then Vt > and x G Ai 

K t (x, x; g) < / + A'(ht/2m)~ D/2 , (64) 



V(M) 

where A' = A'(diam(M),K,V{M)). 

2) Upper Bound of Off-Diagonal Heat Kernel: The following corollary given in [36] 
constrains the off-diagonal elements of the heat kernel from above. 

Assume that for some points x, y G Ai (Ai is any Riemannian manifold) and Vi > 0, 



C C 
K t (x,x;g) < —— and K t (y,y;g) < 



f(t) 



g(t) 



(65) 



where / and g are increasing positive functions on (0, oo) satisfying the regularity condition 
given below. Then, for any C2 > 2 and for all t > 

md 2 (x, y) " 



K t (x,y;g) < 



4A 



exp 



HCot 



y/f{et)g{et) 

where e = e(C2, a), A and a are the constants from the regularity condition below. 
Regularity Condition: There are numbers A > 1 and a > 1 such that 



(66) 



Has) 
/to 



< A 



Hat) 
fit) 



(67) 



for all < s < t. 

By comparing the equations (|64p and (|65p . we realized that the right hand side of (J64 
can be an explicit candidate for the functions / or g in the theorem above. Hence, we 
could have 

+ A'(ht/2m)- 



f(t) = g(t) 



-D/2 



(68) 



V(M) 

by choosing C = 1. It is easy to check that these functions are positive and increasing. We 
can also verify that they satisfy the regularity condition (|67p with A = aPI 2 . Therefore, 
we have obtained the upper bound for the off-diagonal elements of the heat kernel 



K t (x,y;g) < 4A 



1 



+ B(£)(ht/2m) 



-D/2 



exp 



md 2 (x, y) 
hCot 



(69) 



V(M) 

where B(e) = A'e~ D ' 2 . 

3) The Lower Bound of Heat Kernel: 

We have a direct theorem about the lower bound on the heat kernel (cf. theorem 5.6.1 
in [37]). We will just give the statement of the theorem: Let Ai be a complete Riemannian 
manifold with Ric^j > 0. Then, we have 

md 2 (x, y) ' 



K t (x,y;g) > (^M/2m)- D l 2 



exp 



2ht 



for all x, y £ A4 and t > 0. In particular, we find the lower bound to be 

K t (x,x;g) > (4Trht/2m)- D/2 . 



(70) 



(71) 
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4.2 Cartan-Hadamard Manifolds: 



A manifold A4 is called a Cartan-Hadamard manifold [38] if A4 is a geodesically com- 
plete, simply connected, non-compact Riemannian manifold with non-positive sectional 
curvature everywhere. The D-dimensional flat M D and hyperbolic spaces are the best 
known examples of Cartan-Hadamard manifolds. 

1) The Upper Bound of Heat Kernel: In order to give an upper bound for Cartan- 
Hadamard manifolds, we need to give the some definitions and related theorems in the 
literature. 

Isoperimetric Inequalities: Isoperimetric inequalities are the relations between the 
boundary area of regions and their volume. We say that manifold A4 admits the isoperi- 
metric function / if for any precompact open set £1 C M with smooth boundary 

A(dn) > I(v) , (72) 

where A(d£l) is the area of boundary of the region 0, and v = V(£l) is the volume of 
the region. Any Cartan-Hadamard manifold M of dimension D admits the isoperimetric 
function I(v) = kv~ ~ , k > [39] , We have an important theorem [JO] given below: 

Assume that manifold Ai admits a non-negative continuous isoperimetric function 
I(v) such that I(v)/v is non-increasing. Let us define the function f(t) by 

/(*) 



1 = 4 ' dv PM ' (73) 







assuming that the integral does not divergent at t = 0. Then for all x G M, t > and 

e > 0, 

Kt{x ^ 9) -JW^)- (74) 



Moreover, if the function / satisfies in addition the regularity condition (|67p then, for all 
x,y€M,t>0,C2>2 and some e > 0, 

K t ( X ,y;g)< m ex P ^ — j . (75) 

The isoperimetric function for Cartan-Hadamard manifolds given above satisfies all the 
requirements above, that is, it is non- negative continuous function and I(v)/v is non- 
increasing. Substituting the above isoperimetric function for Cartan-Hadamard manifolds 
into (|73p we obtain the function f(t) by a simple integration 

fit) = fa) , (76) 

and it meets all the requirements of this theorem including the regularity condition. Hence, 
the upper bound on the off-diagonal elements of the heat kernel on Cartan-Hadamard 
manifolds is given by 

v( C(e,k) f-md 2 (x,y)\ 

Kt[X ^ 9) ~ {A.Ht/2mr/ 2 6XP { HC 2 t ) ' (77) 

where C(e, k) = T^WjyjDjz an d the physical parameters h, m are introduced for dimen- 
sional reasons. This upper bound estimate of the heat kernel is also valid for mini- 
mal submanifolds, which are submanifolds of M. D whose normal mean curvature vector 



1F> 



H(x) = (Hi(x), H 2 (x), Hf)(x)) vanishes for all x G M, because minimal submanifolds 
admit the same form of isoperimetric function [¥T] with the Cartan-Hadamard manifolds. 

2) The Lower Bound of the Diagonal Heat Kernel: On Cartan-Hadamard manifolds, 
the lower bound of the diagonal elements of the heat kernel are obtained in [38] . The lower 
bound of the diagonal heat kernel is given also for any manifold given as a proposition 
(5.14) in [38]: For any manifold M, for any x € A4 and 5 > 0, there exists c = c x > 
such that 



K t (x,x;g) > 



r exp 



(78) 



(4irht/2m) D / 2 

where a\(M) is the spectral radius corresponding to the Laplacian and it is restricted to 
the following range [42 \ [43] 



(D - l) 2 K 2 max > a x (M) > \{D - l) 2 K 2 min 



(79) 



for a Cartan-Hadamard manifold whose sectional curvature is bounded from above by 
-K 2 ■ . 



5 Pointwise Bounds on Wave function 

There is extensive amount of literature on the exponential decays of the wave functions 
of the Schrddinger operators, which states that L 2 solutions of (—V 2 + V)i/j = Eip obey 
pointwise bounds of the form 

Mr) | < C a e~ ar , (80) 

if the potential energy V is continuous and bounded below and E is in the discrete 
spectrum of —V 2 + V (see [30J for the review of the subject). The proofs given in the 
literature do not include the potentials which require renormalization and they are valid 
only for M. D . We shall prove that it is still possible to get exponential pointwise bounds 
for our problem. 

It is easy to see the upper bound of the wave function ([39]) by applying Cauchy 
Schwartz inequality 



I^fc(a0l < 



< 



< 



a 



a 



N 



Em 

i=l 
- N 

E 



-I) 



dt 



dt tv k 

— e~~ K t (ai,x;g) 
n 



h 



h K t (ai,x;g) 



1/2 



N 







(81) 



where we call the coefficient (which is independent of the coordinates x) in front of the 
integral in (|39p as a and we have used the fact that ^Dl 2 = 1- Thanks to the 

upper bound on the heat kernel given in (169p and (|77p . we show that the wave function is 
pointwise bounded on M. For compact manifolds, the upper bound (|69p of heat kernel 
gives 



N 



\M*)\ < SaA^ 



1 



V{M) 



' md 2 (ai,x) 

v 2 n 2 c 2 



Ki 



I md 2 (a,i, x)v% 

h 2 c 2 



if? 



+ 



B(e) ( md 2 (ai,x)\ 2 4 



h{h/2m) D/2 V ^C: 



Kd 



' md 2 (a,i, x)v\ 



(82) 



where we use the following integral representation of the modified Bessel functions of the 
third kind (or sometimes called Macdonald's functions) 



K v (z) 



2 I 2 



ds e -^V^) s -v-i 



7T 1 

|argz| < -, Re(u) > -- 



(83) 



For D = 2 and D = 3, we can also find the upper bounds on Bessel functions Kq and K\ 
with another useful integral representation 



K u (z) 



Re(z) > 0, Re(v) > -- . (84) 



Using the inequality coshs 
K (x) for x G R + 



ds e~ zcoshs sinh 2,/ s , 

_ e s +e > ^ £ Qr a \\ s y q [ n (jg^l) , we have an upper bound 



/>oo 

K (x) < ds 
Jo 



e 2 e 



By subsequent change of variables £ = e s and r/ = £ — 1, we get 



^o(s) < 



d?/ 



e 2 
?? + 1 



(85) 



(86) 



If we also define a new variable z = xrj, we have 



K {x) < e 2 / dz < 

Jo z + x x 

2 x 

< — e 2 . 



dz e 2 



(87) 



Alternatively, we can find a sharper bound for Kq if we divide the integral above into two 
parts as follows: 



Kq(x) < e 2 

< e 2 



e 2 
dz h 



dz 



e 2 
z + x 



dz 



1 



+ 



1 



Z + X 1 + X 



dz e 2 



Hence, we have the following bound for Kq which shows the logarithmic singularity near 
x = 



Kn(x) < e 2 



< 



In 



x + l\ 2(1 -e" 1 / 2 ) 
1 + x 
x + V 



1 + X 



e 2 + e 2 In 



x 



(89) 



Using sinh s = (- — ^ — ) 2 < for all s > and following the same steps above for K\, 
we find 

i^i(x) <e"i ( - + £ ) • (90) 



s x 2 y 

Substituting (JSS1) and ([90]) into ([82]) for L> = 2, we get 
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N 



\i/> k (x)\ < 8aA^ 



i=l 



l md 2 (ai,x) / md 2 (a ., x)u 2 /h 2 C2 



2V{M)V u 2 h?C 2 



md 2 (ai,x)vl/h 2 C2 



+ 1 



h(h/2m) 



In 



2^md 2 {a u x)v 2 /h 2 C 2 + l 
2^md 2 {a h x)u 2 /h 2 C 2 



l + 2 x lmd 2 {ai,x)vl/h 2 C 2 , 



and for D = 3, we have 



A' 



i=i 



1 



1 / ?^ 2 (a», X) / md 2 (aj , x)t/ 2 /tf > C ,, 

2y(A^)Y ^ 2 C 2 

jB(e)v^ e~ 2 v /md2 ( a »' x )^ 2c ' 2 



+ 1 + 



'md 2 (ai, x)v%/H 2 C 2 
where we have used the explicit exact expression for Ki 



h{h/2mf' 2 2^md 2 {ai,x)/hC2 



(91) 



(92) 



Ki{u) 



(93) 



We can repeat the same steps for Cartan-Hadamard manifolds by using the upper bounds 
of heat kernel given in (I77p and the result is 



N 



\^ h {x)\ < 2aJ2 



t=l 



C(e,k) c - y / md 2 (a ^ x)u 2 /h 2 C2 
h(4irh/2m) 



In 



2Jmd 2 (a il x)v 2 /h 2 C 2 + l 



2Jmd 2 {a i ,x)v 2 /h 2 C 2 



l + 2Jmd 2 {a u x)u 2 /h 2 C 2l 



(94) 



for D = 2 and 



v 



1=1 



^(4vr^/2m) 3 / 2 ( m j2 fax^ftfth) V* 



(95) 



for D = 3. In standard quantum mechanics, the pointwise exponential bounds do not 
take into account singular interactions. Nevertheless, we prove that they are still valid for 
our problem in two and three dimensions. 

In order to understand heuristically why our problem can be considered as a self- 
adjoint extension, which is also suggested by Krein's formula obtained in section [2j it is 
interesting to calculate the expectation value of the free energy for the bound state. The 
result (|39p permits us to write the expectation value as 



{i>k\H \i>k) 



IS 



d 9 x / ~r e h 2l A ^~ u k) K tx{^x;g) / — e » 5J^-( 



A) 



2rr^ 9 



V 2 g K t2 (aj,x;g) 



(96) 



(ipk\H \ipk) 



Using (|17U|) and making an integration by parts to t<i integral, we have 

J2 A U~4)J o J (^K t (a i ,a j ig)e-^A j (-u i 

f°° dti 2 , „ . 

\ -r-e h }_^A*{-ul)K tl {a u x\g) 

+ l^ A ^~ v k) ~T Vke h Kt 2 {aj,x;g) 

„" 1 J 



A' 



(97) 



where we have used the initial condition of heat kernel f)171 1) . Integrating with respect to 
x and using the semigroup property of heat kernel (|172p . we obtain 



(^k\Ho\ip k ) 



N 



5> 



*' 4) 



2 AT 



00 di ft 



K t (ai,aj;g)e ft A,-(-i/ fc ) 



• 7o »,i=i 

~/ it I X S A*(~ l/ i)^i(- l/ I)^i+t2( a i> a i;9)e" 
i,j=l 

By change of variables it = ti + £2 an d i> = ii — we find 

AT 



h vi 



(98) 



(V'jfcl-HolV'fc) 



f°° r\t / t \ tut 



,2 W 







00 du 

7? 



(99) 



One can easily see that the i = j term of the sum for the first term in the parenthesis 

f e-^I^C-i/Dl^Coi.Oi^) (100) 

is divergent due to diagonal short time asymptotics of heat kernel (|18ip for D > 2 and 
one can also show that the off- diagonal terms of this sum in the parenthesis and all terms 
in second sum is convergent by upper bound on the heat kernel given in (|69p and ()77|) . 
Hence we find that the expectation value of the free Hamiltonian is divergent. 



(ip k \Ho\ipk) -> • 



(101) 



If) 



It is a well known fact that point interactions on W D can be considered as a self adjoint 
extension of the free Hamiltonian [2l [10]. We may heuristically think of our problem 
as a kind of self adjoint extension since the wave function ipk( x ) that we have found 
does not belong to the domain of the free Hamiltonian so self adjoint extension of the 
free Hamiltonian extends the domain of it such that the states corresponding to the 
eigenfunctions ipk(%) are included. 

6 Lower Bound of the Ground State Energy E gr 

Although we renormalize the model, we have not completely proven that the energy is 
bounded from below. A well-known theorem in matrix analysis, called Gersgorin theorem 
|44j states that all the eigenvalues a; of a matrix <I> 6 Mn are located in the union of N 
discs 

N 

U{l"~*«l = (102) 

i=i 

where -R£(<&) = X^i=i * s t ne deleted absolute row sums. Since the matrix $ij(E) 
is Hermitian due to the symmetry property of heat kernel, we have u G M. Indeed, all 
eigenvalues co are zero in our problem. If there is a lower bound on energy, that is, ground 
state energy, then we must expect that there would be no solution at all beyond this lower 
bound, say E* = —v*. Then, we want uj = not to be an eigenvalue, thereby none of the 
discs defined above should contain the zero eigenvalue when E < E*. This means that we 
should impose 

N 

| - $u(E)\ = $u(E) > > (103) 

for all i, that is, the principal matrix must be strictly diagonally dominant in order not 
to have a zero eigenvalue. However, before imposing this condition we can simplify the 
problem. We note that 

\$u(E)\ > \$u(E)\^ 

N 

Y;\^ij(E)\ < {N-l)\Q ij (E)r>*, (104) 

so the above condition (I103p is implied by the stronger requirement 

\$u(E)\ min >(N- l)|<D 4J (^)r ax . (105) 

Once we obtain a solution to this inequality, it is satisfied for all E < E* since the diagonal 
part of the principal matrix (126|) is a decreasing function of E and the off-diagonal part 
of it is an increasing function for given etj, Oj and N. This means that there is no solution 
beyond this critical value E*. Hence, the ground state energy must be larger than the 
critical value E*\ 

E gr > E* . (106) 

The basic idea of the proof was given for special manifolds S 2 , H 2 and H 3 in our previous 
work [25]. In the next subsections, we will find the lower bound of the ground state energy 
for more general class of manifolds. 
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6.1 Compact Manifolds 

For compact manifolds, the upper bound for the off-diagonal elements of heat kernel and 
the lower bound for on-diagonal part of the heat kernel has been given in (|69[) and (I71D . 
respectively. Using these bound estimates in (j26j) . we find a lower bound for the principal 
matrix 



1 



(4iTh 2 /2m] 
h(4irh/2m) 



In (u 2 /n 




(107) 



and an upper bound for it 



\*ij(E)\ < { 



4A 




AA 



C 2 V(M)vv d 
2 K i 



+ 



2S(*)*oU/&£ 



(47rn 2 /2m) 




_2_ v_ 

C*2 /-id 



v / 27rC^S(e)/i (J exp 



C 2 K(A*)i/^ d 



+ 



_2__ 
C2 Md 



m 



.3/2 



if D 



if D 



3, 
(108) 

where i ^ j and we have used the monotonic behavior of the functions in &ij(E) so 
that we could maximize the principal matrix in which we defined d = min^,- dij and 
fi = maxj ^ . We also introduced a natural energy scale n d = ^ or simplicity. In 

order to solve the inequality analytically, we must estimate the bounds on the Bessel and 
logarithm functions. We can estimate the lower bound of logarithmic function 



In it > 



u 



1 



for u > , u ^ 1 



Let us first consider two dimensional case. As a result of bounds ([5] 
Bessel functions and the one for logarithmic function given above, we find 



2m 
^2 



1 



1 



and 



\$ij(E)\ < 4A exp 



1 v 



1 



(109) 
90]) given for 

(110) 



V{M) 1 v 2 + ^2C- 2 v^ d 



+ 



2^2C~ 2 B{e) 
(4irh 2 /2m) u/fj, d 



rim 



Since v > [x d , we may have 



$ 4j (£)| < 4A- 



exp 



1 

2C 2 



V2C2~ j {^tf/2m) 
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Therefore, there will be no solution to the eigenvalue equation for values of the ground 
state energy below a critical value v > if the following inequality is satisfied 



2m 
^2 



1 



> AA(N-l)- 



exp 



1 V 

2C 2 Md 



1 



1 + 



1 



V{M)y? d \ V2C: 
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+ 



2^2C^B(e) 
(4vr/i 2 /2m) 



(113) 



Hence, we can solve this and get the lower bound for the ground state energy 

2 



where 



exp 



Egr > -K 



1 n 

2C 2 li d 



-A 



I'd 



+ ^/2C~ 2 w({N - 1)aA 



(114) 



V(M)fi 2 d 



1 + 



+ 



2y^B (e) 
(4vr/i 2 /2m) 



4A(7rh 2 /2m) , (115) 



and W is called Lambert-W function, also called the Omega or the ProductLog function 
l46l and it is defined as the inverse function of xe x . In other words, 



y = xe 



x = W{y) . 



(116) 



As for the three dimensional case, the off-diagonal part of the principal matrix have the 
following upper bound 



$ij(E)\ < 4Aexp 



1 v 



1 



V{M)ti 



1 + 



+ 



^C~ 2 B(e)^ d 
h(47rh 2 /2m) 



(117) 



where we use less sharp upper bound by using v > fi^ in order to solve the inequality. 
Therefore, we conclude that there exists a critical value of bound state energy v > for 
a given N and d such that w^Oso the ground state energy cannot be less than —v 2 



where 



Egr > -v* 



/' 



<C> 



^ + \I^W((N-l)A 2 ) 

H'd 



A, 



exp 



2 y. 
C2 



\f2rtC~2 
h 3 / 2 4A (4irh/2 



V(M)fj& 



1 + 



+ 



V2l^B(e)ix 



m 



3/2 



V2C^ J h?l 2 (47rh/2m) 3/2 



/<■</ 



(118) 



(119) 



6.2 Cartan-Hadamard Manifolds 

Similarly, using the upper and lower bounds of the heat kernel for Cartan-Hadamard 
manifolds, we have obtained the upper and lower bound on the off and on diagonal part 
of the principal matrix, respectively: 



(47rfr 2 /2m) 
2/ttc 
^ h(4iTh/2m 



In 



,3/2 



' V 2 

h 








f. 


\v 2 




h 



]fD = 2 



/r 

n 



(120) 



where £ = h{ai{ ^ )+5) > and 



r 2C( £ ,k) / /Xii 

(47r^ 2 /2m) VV C 2 Md 



I ^ 3 / 2 (4vr^/2 



,3/2 




2 i/ 



if 25 = 2 
if 25 = 3 



(121) 



for i 7^ j. For 25 = 2, we have 



In 



> In 



(122) 



and if we use the same bounds for the Bessel function and the fact that v > fj,d, we obtain 
the lower bound for the ground state energy 



/ 



E gr > -2C 2f i 2 d W 2 



2(N-1)C(e,k) 



\ 



\ 



In 



(123) 



/ 



For three dimensional case, we must do some additional assumption in order to get an 
analytic solution. We will assume that v > p? + fti; and this assumption should be checked 
whether it is consistent or not after we have found the solution. Hence, 



h 1 / 2 




and finally 



where 



Egr > 



-A 



1 

fJ-d 



VpJTh£ + 




W(A 3 (N-1)) 



1 



(124) 

(125) 
(126) 



7 Non-degeneracy and Positivity of the Ground 
State 

The rigorous proof of non-degeneracy and positivity of the ground state in standard 
quantum mechanics is given in |30|. 147] . which does not include the singular potentials. 
Therefore, it is necessary to check whether this is still valid for our problem. The proof 
in our case is based on Perron - Frobenius theorem [44]: It states that If A £ Mn and if 
we suppose that A > (that is, all Aij > 0), then 

(a) p(A) > 0; 

(b) p(A) is an eigenvalue of A; 

(c) There is an x € with x > and Ax = p(A)x\ 

(d) p(A) is an algebraically (and hence geometrically) simple eigenvalue of A; 

(e) < p(A) for every eigenvalue u / p(A), that is, p(A) is the unique eigenvalue 
of maximum modulus. Here p(A) = max{|w| : uj is an eigenvalue of ^4} and called 
spectral radius. A simple proof of this theorem for the positive symmetric matrices is 
given in |48j . 



5>3 



Since the principal matrix (|26p is not a positive matrix, we cannot directly apply 
Perron-Frobenius theorem. Nevertheless, we can make the principal matrix strictly pos- 
itive by subtracting the maximum of the diagonal part of it corresponding to the lower 
bound of energy E = E*, which is found in section [6l and reversing the overall sign: 

$'(£) = -[$(£) -(l + e)l$£ ax (£*)]>0, £>0, (127) 

where 1 is a N x N identity matrix. Hence, considering the transformed principal matrix 
<3?' in the light of this theorem, we conclude that there exist a strictly positive eigenvector 
which corresponds to the unique eigenvalue of maximum modulus. 

N 

Y,^( E )ME) = P{&)ME). (128) 

3=1 

Here it must be noticed that has the same eigenvector with so it guarantees that 
there exist a strictly positive eigenvector for the principal matrix <J>. Using the eigenvalue 
problem, we find 

N 

J2^ij(E)A j (E) = uj{E)A t (E) , (129) 
i=i 

where p($') = — u mm (E) + (1 + e)&u(E*). For a given E = Ek or v = z^, there is a unique 
corresponding u> mm (E) and since we are looking for the zeros of the eigenvalues co(E) = 0, 
this minimum flows to zero at v = z/ nax = . This means that the positive eigenvector A± 
corresponds to the ground state energy so we prove that the ground state energy is unique 
and the associated eigenvector Ai is strictly positive. Due to the positivity property of 
heat kernel, it is easy to see that the ground state wave function is strictly positive from 
the equation (|39|) . 




(130) 



Hence, we prove that despite the singular character of the interaction, the ground state 
is still non-degenerate and unique. 

8 Renormalization Group Equations 

Renormalization group equations of our problem (for N = 1 case) in flat spaces has been 
already given in the literature [6| ll5lfT8] . Here, we will show that this can be also derived 
explicitly for our problem. 

8.1 Two Dimensional Case 

In this section, we shall choose the natural units h = 2m = 1 for simplicity. It is useful to 
work with the dimensionless coupling constant. One possible way for the renormalization 



OA 



scheme in order to determine how the coupling constant changes with the energy scale 
is to define the following renormalized coupling constant A^(Mj) in terms of the bare 
coupling constant Aj(e) 



Af(M0 Xi(e) 



-Mft 



dt- 



4vrt 



(131) 



where Mj is the renormalization scale (it is of dimension [i?] 1 / 2 ). Then, the renormalized 
principal matrix in terms of renormalized coupling constant in natural units is given 



mm 



i 



dt K t (a i ,a i ;g)e t 



-M?t 



) 

dt Kt(ai,aj;g)e tE 



4vrt 



if i = j 
if i£ 3 , 



(132) 



and the bound state energy is determined from the condition det^ R (E) = and it 
determines the relation between A^(Mj) and Mj. Explicit dependence on M, cancels the 
implicit dependence on Mj through \f-(Mj). Physics is determined by the value of \f(Mi) 
at an arbitrary value of the renormalization point Mj. However, this is not a proper way 
to look at our problem since we have to deal with several renormalized coupling constants 
with the same kind of interaction, which essentially differs from one another with arbitrary 
constants. These constants can be determined by deciding the excited energy levels. We 
shall instead prefer one renormalized coupling constant by redefining the meaning of 
the renormalized coupling constant without altering physics. This could be done in the 
following way: As an external input, we decide about the relative strengths of individual 
delta interactions and do not use the ground state energy to fix the flow. We know that 
—fj% is the bound state energy of the individual i th Dirac delta center so it corresponds 
to the solution yu 2 ) = 0. Without loss of generality, let us assume that $£(-^1) = 
and this allows us to choose the renormalized coupling constant 



A«(M) Ai(e) 



-m'H 



dt- 



4vrt 



(133) 



at some scale M. Once the renormalized coupling constant is fixed under this condition, 
we must also satisfy fJ-f) = for i 7^ 1 with this choice at the same scale M. This 

is always possible if we add a constant term to the definition of renormalized coupling 
constant. Let us consider i = 2 case 



1 



. + lim 

X R (M) e^O 



00 e ~ M2t 



dt 



4vrt 



dt K t (a 2 ,a 2 ;g)e ^ 



dt 



K t (a 1 ,a 1 ;g)e ^ - K t (a 2 , a 2 ; g)e 



(134) 



where we have used (T133D and $f x (-/x 2 ) = 0. This means that there always exist a 
constant Sj depending only on /Xj with Si = and Sj 7^ for i 7^ 1 such that the 
condition $^(— /U 2 ) = can be satisfied. Hence, the renormalized coupling constant 
becomes 



1 



1 



X R (M) Aj(e) 



-M 2 t 



dt 



Ant 



(135) 



and the choice of Ej s refer to the relative strengths of delta interactions in this renormal- 
ization scheme. As a result, the renormalized principal matrix is 



I , . / „ , s tE e 



dt K t (ai,aj;g)e tE if j. 

o 



(136) 



The renormalization condition is given by 



dM 

or 



d$^(M, X R (M),E;g) 
M l]K ' "\ ' ' = , (137) 



d , d 



R 



where 



Mr) = m ^m (139) 

is called (3 function and the equation (|138p is the renormalization group (RG) equation. 
In [15], the renormalization condition (|137p corresponding to the problem in flat space 
has been written in terms of the T-matrix. Using (|136p in (|138p , we can find /3 function 
exactly 

/3(A R ) = < . (140) 

This result is the same as the one in flat spaces given in the literature [UJ [T4"l [To"! [15] so our 
problem is asymptotically free. From the explicit expression of the renormalized principal 
matrix, one can easily see the scaling property of it under a change of energy and metric 
scale 7 using the scaling property of the heat kernel (|179p 

X R (M),^ 2 E; j~ 2 g) = $g(7 _1 M, X R (M ), E; g) . (141) 

It is important to notice that we need to scale the metric as well and the idea of the 
metric scaling in deriving the renormalization group equation was motivated by [49] in 
the context of renormalization group in quantum field theory on curved spaces. Hence we 
have 

7-^ [$g(M, X R (M), 7 2 £; T 2 g) = *g(7 _1 M, X R (M), E; g)] . (142) 
This leads to the renormalization group equation for $g(M, X R (M), 7 2 £; 7 " 2 5 ) 

1 ^f j (M,X R (M), 1 2 E ]1 - 2 g) + M^^(M,X R (M), 1 2 E ]1 - 2 g) = 0, (143) 

or 

7^"/3(A R )^ ^(M,X R (M), 1 2 E; 1 - 2 g) = 0. (144) 
If we postulate the following functional form for the principal matrix 

<P? j (M,X R (M), 1 2 E;r 2 g) = f(j)^(M,X R (jM),E;g) , (145) 
and substitute into (|144p we obtain an ordinary differential equation for the function / 

7^ = . (146) 
d7 

9fi 



This gives the solution 7(7) = 1 using the initial condition at 7 = 1. Therefore, we get 

<^(M,A fi (M),7 2 £; 7 - 2 <7) = ^(M, X R ( 7 M), E; g) , (147) 
which means that there is no anamolous scaling. After integrating 

between M = M to M = 7M we can find the flow equation for the coupling constant 

A «^ M > " ,4. "frlni ' (149) 

1 + 2^A jR (M)ln7 

One can explicitly check the relation (|147p if the coupling constant evolves according to 
(11491) , First, we add and subtract a term in the time integral: 

11 f°° ( e~ M2f \ 

^(M,X R (lM),E; 9 ) = — + -l^-j Q dt ^K t (a h <H; g)e* E - — J - 

1 1 , f°° '( , , ^ e~ MH e -^V 2 t e - M ^~ H \ 



rco I e -M "t 



A K (M) 2tt ' J I v ' 4vrt Ant 4irt J 

1 f°° { e -M 2 7 - 2 t\ 

dt K t (a hai -g)e tE — — - Si (150) 



and then using the scaling property of heat kernel (|179p . we get 

dt [ 7 _2 i^ 7 -2 t (ai,aj;7 _2 5)e* 



1 f 00 / 2 e -A:f2s \ 

A^(M) " 7o ^ [X'faW 2 ^ E -^)- Ei 
$%(M,\ R (M),^E; 7 - 2 g). (151) 



Off diagonal term can be directly checked using just the scaling property of heat kernel 
(|179p . In another way of thinking, one can find how the coupling constant evolves, given 
(|149|) from the scaling relation (|147p . 

8.2 Three Dimensional Case 

Since it is convenient to work with the dimensionless coupling constant, we define a 
dimensionless coupling constant in three dimensions 

Xr(M) = MX R (M) . (152) 

Then, by similar arguments developed for two dimensions, the renormalized principal 
matrix in the natural units is 



\ R {M) Jo V (4vrt)3/ 2 y (153) 

- \ dt K t (ai,aj;g)e tE iti^j- 
Jo 



9,7 



Renormalization condition (|138[) in this case leads to the following (3 function 



(154) 



which is in agreement with the result for flat space [15]. From the explicit expression 
of the renormalized principal matrix for three dimensions, one can easily see the scaling 
property of it under a change of scale 7 using the scaling property of the heat kernel (|179p 



^ j (M,X R (M), 1 2 E; 1 - 2 g)= 7 ^( 7 - 1 M,X R (M),E;g) , 



so we have 



(155) 



7 



d7 



$g(M, A fl (M), 7 2 £; 7" 2 5) = 7^(7 _1 M, X R (M),E; g) 



(156) 



This leads to the renormalization group equation for $^(M, A\r(M), 7 2 -E; 7 2 <?) 



7— -1 + M 
d7 



$? i (M,\ R (M), 1 'E; 1 -*g) = 0, 



or 



<9M_ 

d 7 <9A fl 



^(M,A /? (M),7^;7- 2 5 ) = 0. 



(157) 



(158) 



If we again postulate (|145p and substitute it into (|158p . we obtain an ordinary differential 
equation for the function / 

d /(7) , ,,- a s 
7 — = f . (159) 

The solution is 7(7) = 7 by using the initial condition at 7 = 1. Therefore, we have 

^(M,X R (M), 1 2 E- 1 - 2 g)= 7 ^(M,\ R ( 1 M),E;g) . (160) 



This means that there is also no anamolous scaling in three dimensions. After integrating 

(161) 



= M 9 -^ 1 = UM) - ^Ai(M) , 



between M = M to M = we can find similarly the flow equation for the coupling 
constant 

7 A fl (M) 



Xr^M) 



(162) 



l--LA R (M)(l- 7 ) 

One can similarly check the relation (|160p if the coupling constant evolves according to 
(|162p . In this case, we have 
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-M 2 7 - 2 i\ 

M I KJn,. n.- n\p m - — 

MM) Jo 

and then using the scaling property of heat kernel (|179p . we get 



roc / g _M 2 7 - 2 A 
7 Jo ^ \ Kt ^ a ^ 9 ^ tE - (4^)3/2 J 



M f°° I e -M 2 7~ 2 i 

7/ dt [j-ZK^tia^ain-ig^ 



A fl (M) Jo V (4vrt) 3 / 2 
M 



j ds ( K s (a i ,a i \'y 2 g)e 



-M 2 s 



(163) 



= <^M,A R (M), 7 2 £; 7 - 2 5 ) • (164) 

One can similarly find how the coupling constant evolves, given in (|162p from the scaling 
relation (TTdOI) . 



9 Conclusion 

In this paper, we studied the bound state problem for several Dirac delta interactions in 
various two and three dimensional Riemannian manifolds. Our renormalization method 
is basically inspired from [23] developed for many body version of this problem on flat 
spaces. This method allows us to renormalize the model non-perturbatively. It has been 
also shown that the heat kernel plays a key role in the renormalization procedure and 
help us to find lower bounds on the ground state energy due to the sharp upper bound 
estimates on it for several classes of manifolds. We proved that many well known theorems 
given in standard quantum mechanics are still valid, such as pointwise bounds on the wave 
functions, non-degeneracy and uniqueness of the ground state although we have singular 
interactions. Renormalization procedure does not change these well-known results in 
standard quantum mechanics. Finally, we studied the renromalization group equations 
and the /3 function is exactly calculated. 
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11 Appendix A: Heat Kernel on Riemannian Man- 
ifolds 

Let (Ai,g) be a compact connected Riemannian manifold, then there exists a com- 
plete orthonormal system of C°° eigenfunctions {//}^ in L 2 (Ai,d^x) and the spectrum 
a(M.,g) = {ai} = {0 = uq < o\ < 02 < ■ ■ ■ }, with 07 tending to infinity as / — > 00 
and each eigenvalue has finite multiplicity: —V 2 g fi{x;g) = (Jifi{x;g) . The multiplicity of 
the first eigenvalue <jq = is one and corresponding eigenfunction is constant and given 
by fo{ x '-,g) = 1/ \/V(A4), where V(Ai) is the manifold (A4,g). This theorem is also 
called Hodge theorem for functions or spectral theorem [43 \ [50] . This theorem is valid for 



'X) 



Neumann, Dirichlet and mixed eigenvalue problems except for ctq > 0, provided that the 
appropriate boundary condition is imposed. 

The operator — V 2 is formally self-adjoint with respect to the L 2 (Ai, dfx), that is, 
the inner product is defined as 

(fo,ih)= I dfx^l{x)^ 2 {x) • (165) 

Jm 

Since it will be necessary for our purposes to know how the eigenfunctions change under 
the scaling transformations in the metric, we will use a notation which specifies the metric 
structure of the functions, such as fi(x; g). The spectral theorem provides us with all the 
tools of Fourier analysis, so that we can expand any function ifj(x) G L 2 (A4, df x) in terms 
of the complete orthonormal eigenfunctions fi(x;g) 

ip(x) = (^)Ji(x; g)) fi(x; g)=Y, c i M*; a) < ( 166 ) 

l>0 l>0 

Orthogonality and completeness of the eigenfunctions on compact manifolds are 

Ski = dfx fk(x;g)fi(x;g) , 
Jm 

5g(x,o,i) = ^2fi{x;g)f?{ai;g), (167) 
z>o 

where Cj's are expansion coefficients. 5®(x,ai) is the D - dimensional normalized delta 
function at point a, G M, 

[ tfx6f(x,ai) = l. (168) 

JM 

Note that extra labels in the eigenfunction expansion must be taken into account if the 
problem admits degeneracy. 

Moreover, one may extend heuristically the problem on to some noncompact manifolds 
in such a way that spectral theorem is still applied. The relations such as completeness 
and orthogonality relations are defined by an appropriate generalization of the measures 
to the continuous distributions in the sense of |51| . Since the spectrum does not have to 
be discrete for noncompact manifolds, we may have in general 

i>(x) = j Mi) mfi&g) , 

Sm = dfx fk(x;g)fi(x;g) , 
Jm 

S°(x, ai ) = Jd/iQ) Mx;g)ft( ai ;g), (169) 

where dfj,(l) is the spectral measure and it includes continuous spectrum as well as point 
spectrum. This should be taken with a grain of salt and one must not forget that we may 
not have a rigorous proof of the spectral theorem for some special noncompact manifolds. 
From now on, we assume that we are dealing with manifolds which do not have such 
pathologies, that is, the spectral theorem is applicable. 

Although the notion of heat kernel can be defined on any Riemannian manifold, the 
explicit formulas only exist for some special class of manifolds, for example, Euclidean 
spaces R D [38] and hyperbolic spaces M D [52] , We will list some of the well known 
properties of the heat kernel on any (A4,g), and give the short time asymptotic expansion 



of the heat kernel. Also we derive the scaling property of the heat kernel. They are all 
needed in our calculations. 

1) Heat Equation: It satisfies heat equation since it is a fundamental solution to it by 
definition. 

dK t (x,y;g) fr 2 2 
h j t —V g K t (x,y;g) = , (170) 

where Vjj acts on the x variable. 

2) Initial condition: It solves the Cauchy problem. 

hm i K t (.,y;g) = 5 g {.,y) . (171) 

3) Semi-group Property: 

/ dgZ K tl (x, z; g)K t2 (z, y; g) = K tl+t2 (x, y; g) . (172) 
Jm 

4) Symmetry Property: 

K t (x, y; g) = K t (y, x; g) . (173) 

5) Positivity Property: 

K t (x, y; g) > for all t > . (174) 
If M is compact, then we have 

K t (x,y-g)=Y,e~^ a 'ft(x;g)fi(y;g) , (175) 

l>0 

which converges uniformly on Ai x Ai . The analog of this sum on noncompact manifolds 
can be given as 

K t (x,y;g) = [ d/^e-^Vf 3) • (176) 



We have also eigenfunction expansion of the creation and annihilation operators 

4(*)=2>J(0/ffo(7) (177) 



l>0 

or in noncompact manifolds 



<to= / dfi(l)4(l)ft(x;g) (178) 



When the manifold Ai is a complete Riemannian manifold with Ricci curvature bounded 
from below then the heat kernel satisfies the stochastic completeness property \38\ 153] : 

/ dfxK t (x,y;g) = l 
Jm 

On a compact manifold stochastic completeness is always satisfied [5l]. Using the prop- 
erties of heat kernel and stochastic completeness, one can derive the scaling property of 
the heat kernel 

K t (x, y; g) = a D K a 2 t {x, y; a 2 g) . (179) 
Free resolvent can be written in terms of the heat kernel 

1 f°° 

R Q (x, y\z) = (x\ (H - zy 1 \y) = -J dt e^K t (x, y; g) , (180) 

31 



where Re (z) < 0. It can be defined for complex values z by analytic continuation. We 
have short time asymptotics of the diagonal heat kernel for any manifold [55] 



1 

^ ( ™ } ~ (4^/2^ g Ufc(a ' J;)(fa/2m)fc/2 ' (181) 

for every x £ M. Here the functions Uk(x,x) are scalar polynomials in curvature tensor 
of the manifold and its covariant derivatives at the point x. When there is no boundary, 
the odd terms in the expansion, i.e, k = 1, 3, 5, ... vanishes |56j . In this paper, we always 
assume that the manifolds have no boundary. Indeed, we have also short time asymptotic 
of the heat kernel for any x and y with the several assumptions about the structure of 
the set of geodesies which join the points x and y |57j . It is shown that (see Theorem 2.1 
and 2.2 in [57]) for all y sufficiently close to x (so that x and y can be joined by a unique 
shortest geodesic j XtV along which x and y are non-conjugate) then, 



^'»~ (iap^ ' v ^'' ) ' (182) 

where d(x,y) is the geodesic distance between x and y and ^ 7 (x,y) characterizes the 
divergence of the geodesic flow near 7, that is, if we emit a beam of geodesies from x 
along 7 in the solid angle d<p illuminating a hypersurface of area dS at y orthogonal to 7, 
then fy(x,y) = dS/dip. The function *$>(x,y) can also be written in terms of the Jacobi 
fields orthogonal to the geodesic Jx,y(s), where < s < d(x,y). If the number of shortest 
geodesies joining x and y is greater than 1, or, if x and y are conjugate along some of 
them, then the result takes the following form (up to a bounded factor) 

/ (D + k) md 2 (x,y)\ 

K t (x,y;g)=0[(M/2m) ~e ^ (183) 



where the index k = k(x,y) depends on the character of the degeneracy of the geodesic 
flow between x and y. Since all our calculations essentially give the same physical result 
for all cases and subcases [57], we will just consider a generic case (Case 3.1 in [57]): The 
set £l x ,y consists of finite number of geodesies 71, ...,7 m and x and y are non-conjugate 
along each of them. In this case, for each 7, we can define ^?i(x,y) by considering the 
Jacobi fields along 7j. Then, we have (Theorem 3.1 in [57]): 



Kt{X ^ 9) ~ (4^/2m)^ d(g " 1)/2(a; - y) E^ 1/2 (^y) ■ (184) 

12 Appendix B: Existence of Hamiltonian 

Let A be a subset of the complex plane. A family J(E), E S A of bounded linear operators 
on the Hilbert space "H under consideration, which satisfies the resolvent identity 

J(E 1 ) - J(E 2 ) = (E 1 - E 2 )J(E 1 )J(E 2 ) (185) 

for Ei , Ei G A is called a pseudo resolvent on A [58]. The following corollary (Corollary 
9.5 in [SH]) gives the condition for which there exists a densely defined closed linear 
operator A such that J(E) is the resolvent family of A: Let A be a unbounded subset 
of C and J(E) be a pseudo resolvent on A. If there is a sequence E n G A such that 
\E n \ — > 00 as n — > 00 and 

lim -E n J{E n )x = x , (186) 



39 



for all x £ H, then J(E) is the resolvent of a unique densely defined closed operator A. 
In order to show the resolvent kernel that we have found in (|2ip corresponds to a unique 
densely defined closed operator H, we need to prove that it satisfies the resolvent identity, 
i.e, 

R(x,y\E 1 )-R(x,y\E 2 ) = (E 1 -E 2 ) [ dfz R(x, z\Ex)R(z, y\E 2 ) . (187) 

JM 

Substituting (f2Tj) into ()187p . we obtain 

N 

R (x,y\E l )-R (x,y\E 2 )+ ^ Rq{x, a^Ex)^ (Ex)^^ y\Ei) 



N 



(E x - E 2 ) [ d D z 
J M 



- R (x,a i \E 2 )^T. 1 (E 2 )R (a j ,y\E 2 ) 



Ro(x,z\Ex)Ro(z,y\E 2 ) 

N 

+ Ro(x,z\E 1 )R (z,a l \E 2 )^- j 1 (E 2 )R (a j ,y\E 2 ) 

N 



N N 
i,j=l k,l=l 



xR (z,a k \E 2 )$- l (E 2 )R (a h y\E 2 ) 



(188) 



Using the formula ()18Qj) . it is easy to see that the free resolvent satisfies the resolvent 
identity 

(Ex~E 2 )f d^zR (x,z\E 1 )R (^y\ E 2) 
Jm 

= (Ex -Ei)J m dfz ^K tl (x, z; g)<* E ^ ^K t2 (z, y; g)e^ 
= (Ex - E 2 )j^j™ ^^K h+t2 (x, y] g)e^e^ 



(Ex - E 2 ) r°°du 
2 Jo n 



u dv 



K U (X, y, g ) e (n+v)E 1 /2h e (u-v)E 2 /2h 



= ^K u (x,y;g) (e uE ^l h - e uE >' h ) = R (x,y\Ex) - R (x,y\E 2 ) , (189) 

where we have used the semigroup property of heat kernel (|172p and made the change of 
variables u = t\ + t 2 , v = tx — t 2 . Then, the equation (|188p becomes 



v 



]T Ro(x, a i \Ex)^ l (Ex)Ro(a v y\E x ) 



N 



Y Ro(x,a i \E 2 )^ 1 (E 2 )R (a j ,y\E 2 ) 



*>J=1 



(Ex - E 2 ) [ d^; 
Jm 



N 



]T R (x,z\E 1 )R (z,a l \E 2 )^ 1 (E 2 )R (a :j ,y\E 2 ) 
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N 



+ ]T R (x,a i \Ex)^ 1 (Ex)R {a j ^\E 1 )R Q (z, y \E 2 ) 



N N 
i,j=l k,l=l 



xR (z,a k \E 2 )^- 1 (E 2 )R (a l ,y\E 2 ) 



(190) 



If we add and subtract the terms Yli j=x Ro( x , c L i\Ex)&ij 1 (Ex)Ro(a,j, y\E 2 ), 
Y^^j = iRo(x,ai\Ei)^j 1 (E 2 )Ro(aj,y\E 2 ) on the left hand side of (|19Up . and rearrange all 
the terms, we obtain 

N 

RofaoilEjQ^E!) [R ( aj ,y\Ei) - R^a^y^)] 

N 

+ Ro(*,<h\Ei) l^HEi) - <b-\E 2 )\ Ro( aj ,y\E 2 ) 
N 

+ ^ [Ro(x,<k\Ei) - Ro(x,a i \E 2 )]^T. 1 (E 2 )R (a j ,y\E 2 ) 

- r 

= (E 1 - E 2 ) Ro{x,<H\Ei)*J±{E{) / dfz Ro{aj,z\Ei)Mz,y\E2) 
i,j=l J M 

N 

+ Ro(*,<H\Ei) [^(Ei) ~ ^-\E 2 )\ R ( aj ,y\E 2 ) 

- r 

+{E l -E 2 )Y ^zR Q {x,z\E l )R Q {z,a l \E 2 )^- 1 {E 2 )R {a j ,y\E 2 ) , (191) 



where we have used the result (|189p in the first and third terms. The second term can be 
written as 

N N 

£ 2 ^a^i)^ 1 ^^ . (192) 

i,j=l k,l=l 

It is important to notice that difference of the principal matrix equals to the difference in 
free resolvent kernel, that is, 



<f>ij(E 2 ) - <$>ij{Ei) = R {a i ,a j \E x ) - R (a t , aj \E 2 ) 

= {Ex - E 2 ) [ dfz Rofazl^Ro&ajlEz) , 
Jm 



(193) 



for any i and j. After substituting (|193|) into (|192p . the equation (|19ip becomes 

N 

(Ek-Ez) Y R^a^Ex^T^Ex) / d D g z R (a 3 , z\Ex)R (z, y\E 2 ) 



+(Ex-E 2 ) Y E iJoCs.Oil^)^ 1 ^! 

ij'=l fc,Z=l 



: / dfz^ (ai,2|^i)^o(^,afc|^2)^ 1 (^2) J Ro(ai,2/|^2) 
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N 

+(E 1 -E 2 )Y / / dfz R (x,z\E 1 )R (z,a i \E 2 )^r. 1 (E 2 )R (a v y\E 2 ) , 
u=i Jm 



(194) 



This is exactly equal to (|190p so resolvent identity is satisfied. We must now imply the 
following condition in I? norm 



\E n R(E n )f + f\\^0, 



(195) 



as n — > oo and / belongs to the Hilbert space H under consideration and the norm is 
taken with respect to %. Let us choose the sequence E n = —uEq since the resolvent is well 
defined in this resolvent set, in which we have no spectrum below the absolute value of 
the bound Eq that we have found for the ground state energy. Without loss of generality, 
we can set Eq = c\EJ, where c > 2. Then, we have 



\\nE o R(-nE )f-f\\^0, 
as n — > oo. Using (|2ip and separating the free part, we get 



(196) 



\\nE R(-nE )f - f\\ < \\nE R (-nE )f - f\\ 

+ nE \\R (-nE )^ 1 (-nE )R (-nE )f\\ . (197) 

Since it is well known that the first part of the sum converges to zero as n — > oo, that 
is, free resolvent defines a densely defined closed operator (Laplacian), we are going to 
investigate only the second term 

nE \\Ro(-nEo)^-H-nE )R (-nE )f\\ 

N 



< nE { 



^ / dfx R (ai,x\ - nE )R (x,ai\ - nE ) 
j,k,i=i jM 

[ d^yRo( aj ,y\ -nE )R (y,a k \ - nE^r. 1 (-uEq)^ 1 (-nE ) 
J M 



1/2 



,(198) 



where we have used the fact that the Hilbert space norm of an operator is smaller than 
its Hilbert- Schmidt norm: \\Af\\ < Tt 1 / 2 {A^ A) with A = Roi-nEo^^i-nEo) 
Ro(-nEo). By using (|180j) and the change of variables u = t\ + 1 2 , v = t\ — t 2 , we get 



r roc rc 

/ dfx R (a i ,x\-nE )R (x,a l \-nE ) = \ 

Jm Jo Jo 

00 dt *bq 
t K t (ai,ai;g)e 7 



dhdt 2 _ n (ti+t 2 )E 
K tl+t2 (ai,ai\g)e » 



h 2 



h 2 



(199) 



Let us first consider the diagonal case i = I and k = j in the above. Then, the equation 
(fT98|) becomes 



nE 



N 

/ dfx R (ai,x\ - nE )Ro(x,ai\ - nE ) 



x / dfy Ro{ aj ,y\ - nE Q )R Q {y, aj \ - nEo)]^ 1 (-nE )\\^jl(-nE )\ 



M 



1/2 



.(200) 



In this case, the upper bound of the equation (|199p can be found from the upper bound 
of the heat kernel for compact manifolds 



/ dfx R ((ii,x\ - nE )R (x,ai\ - nE ) 
Jm 
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< 
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+ 



AAB(e) 



nEr 



V{M)E 2 n 2 h 2 {h/2m) D / 2 \ H 
and for Cartan-Hadamard manifolds 

/ d^x Ro(cii, x\ — uEq)Ro(x, ai\ — nE$ 
Jm 

C(£,k) 



r 2 



D 
~2 



(201) 



< 



nEr 



h 2 (4irh/2m) D / 2 \ H 



r 2 



D 
~2 



(202) 



In order to give the upper bound for the inverse principal matrix, we decompose the 
principal matrix into two positive matrices 



$ = D - K 



(203) 



where D and K stand for the diagonal and off diagonal parts of the principal matrix. 
Then, it is easy to see $ = D(l — D^ 1 !^). The principal matrix is invertible if and only 
if (1 — Z? -1 /^), and (1 — D~ 1 K) has an inverse if the norm ||D -1 .ftT|| < 1. Then, we write 
the inverse of $ as a geometric series 



= (1 - D- 1 Ky l D- 1 

= (l + iD^^ + iD^K) 2 + ...) D- 1 , 



(204) 



where we must have ||-D -1 -fT|| < 1. Since we are not concerned with the sharp bounds 
on <i> _1 for this problem, we can choose ||-D -1 -ff|| < 1/2 by adjusting the uEq sufficiently 
large without loss of generality and get 



l^" 1 ! < 2\D- 



(205) 



The lower bound of the diagonal principal matrix for compact (1107|) and Cartan-Hadamard 
manifolds (|120p gives the upper bound of the inverse principal matrix. Hence, we find 



-nEn)\ < < 



(4Trh 2 /2m) In" 1 (nE /fi 2 ) 
h{A^h/2mf /2 



2^ 



nE n 



if D 



V- I if D 



(206) 



for compact manifolds and 



m\-nE )\ < < 



{ (4 nH 2 /2m) ^ + t 



h{A-Kh/2m 
20Fc 



,3/2 



nEo 

h 



if D 



(207) 



for Cartan-Hadamard manifolds. If we substitute the results (|20ip and (|202p into (|200p 
for D = 2, and take the limit n — > oo, the result goes to zero. Since the norm is always 
positive, we prove 

\\nE R(-nE )f- f\\^0 (208) 

as n — > oo. It is almost evident that the off diagonal terms in the sum also vanishes in 
this limit primarily because these terms are exponentially damped e _v/ ™ due to the upper 
bound of the heat kernel. Unfortunately, the proof for D = 3 is more subtle and the 



3fi 



volume growth conditions of the manifolds are very delicate in the analysis so we can not 
prove it for three dimensional case by the same approach. Instead we follow the following 
approach: In three dimensions, estimating the Hilbert space norm by Hilbert-Schmidt 
norm does not lead to zero. Instead we will first show that the last term 



\E„ 



[ d 3 g x Mx,ai\En)^{ E n) I d 3 g zR ( aj ,z\E n )f*( 

J M - „■ u i—i JM 



i,j,k,l=l 



R (x,a k \E n )<l> kl l (E n ) / d 3 y R (ai,y\E n )f(y) 
JM 



1/2 



(209) 



goes to zero as E n — > — oo for any / 6 L 2 {M). Prom our previous argument, we know 
that the inverse of the principal matrix <3? satisfies: 



max,- 



< 



A, 



E I 1 / 2 



(210) 



where we define all the constant terms coming from the bounds of the heat kernel as A2 
(exact form of the constants is not important here) and ignore the term in the denominator 
for large values of n for simplicity, which can be read from f|206[) and (|207p . We shall use 
the notation for the constants coming from the bounds of the heat kernel combined with 
the other constants factors as A\, A2, A3, . . . for simplicity. Moreover, we can combine 
the two resolvents with the common variable x, and as a result, we can express this 
combination as 

1/2 



dt t 



\E n \t 

ft K t (ai,a k ;g) 



(211) 



and pull it out of the square root. Using similar arguments as before, we can show that 
this term in three dimensions, for both Cartan-Hadamard type manifolds and Compact 
manifolds (bounded Ricci), including the identical beginning and end points, is smaller 
than 

A 3 



LEU 1 /* 
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where ^3 can easily be read from the upper bound of the heat kernel. Hence we end up 
with the fact that the expression (|209p is smaller than, 



N 



NA A \E n \y* J2 

3,1=1 



d 3 y R Q { aj ,z\E n )\f{z)\ / d;!z #0(2/^)1/(2/)! 



M 



M 



1/2 



(213) 



Hence we should show that the term 



M 



d 3 g yRo( aj ,y\E n )\f(y)\ 



(214) 



decays faster than l^l -1 / 4 . 

To do this, we will pick any one of the centers and choose Riemannian normal coordi- 
nates around it, we assume that the injectivity radius of the manifold is 5 > 0. But first 
we reexpress this term in terms of the heat kernel and use some bounds, 



R (aj,y\E n ) 



°° dt \En\t 

—e n K t (aj,y;g) 



■AT 



for Cartan-Hadamard manifolds, and for compact manifolds we have a similar term with 
an inverse volume term, , added. In the case of compact manifolds volume contri- 
bution term goes to zero faster than ji^l -1 / 4 as can be checked easily, so it causes no 
problems. In both cases we will concentrate on the least convergent part. If we evaluate 
the integral over t now we find, 



r r n / m d 2 (aj,y)\E„\ \f( \\ 

/ d 3 g y R ( aj ,y\E n )\f(y)\ < A G / d 3 g y e V ^ JJU/2L . ( 2 15) 
We divide the right hand side as, 



_ 2 ^^ ]m r -2 ^ ( ^ En] \m\ 



J 

here the last term is smaller than 



d 6 e v ^c 2 jp^i + d t e V ^ 2 jp^i (216) 

(a 3 ) d(aj,y) JM\B s (a 3 ) d{aj,y) 



^ 2 \E n \ 



~h T C 2 f _ / md 2 (aj,y)\En 

I 3 

M\B 6 (a 3 ) 



d^e V \f(y)\ 



e 



'T,S 2 \E n \ 



h 2 C 2 f _ I md 2 (aj,y)\E n \ 

I M 



e 

< — 



md~\En\ , -i In 

md 2 (a,,y)\E n \ 1 i / Z 



h 2 C 2 



dlye 2 V 

M 

By a theorem of Gaffney [59J . for a stochastically complete manifold, for any a > (with 
the inverse length dimension) and any point a on the manifold, the integrals satisfy, 

d 3 e -ad(a,y) < ^ 

IM 

This establishes that the last term decays faster than |S n | _1//4 , so we should look at the 
first part. As has been said we go to the Riemann normal coordinates of the geodesic ball 
of radius <5 and write the integral in terms of Gaussian spherical representation: 



dn drr 2 J(r,6)e V h ' c 2 LLL2 — IL _ (218) 



o ' 

Let us recall that in Gaussian spherical coordinates, the integral of a function / on an 
D-dimensional Riemannian manifold Ai becomes, 

r rpn 

d£xf{x)= dn drr D ~ 1 f{r,6)J(r,9). (219) 

M JS - 1 JO 

Here f2 denotes the direction in the tangent space around a point that we choose, and po, 
refers to distance to the cut locus of the point in the direction Q. We will now divide the 
integral over r to two parts, 



Js 2 Jo r J§ 2 J A 



r 

(220) 



: ; !S 



where < A < 5. Let us now introduce the following function: 



sn fc (r) = < 



sin 

s/k 

r 

sinh(v / — kr) 



if k > 
if k = 
if fc > 



(221) 



This function is very useful for the Bishop-Gunther volume comparison theorems. We 
assume that M has Ricci tensor bounded from below by k\, i.e. Ric(.,.) > kig(.,.) 
and sectional curvature K bounded from above by k<i- Then, the Jacobian factor of the 
Gaussian spherical coordinates satisfies an inequality as follows 



sn? (r) snl (r) 



(222) 



In the second integral, we use 

rS 



dr r e 



< 



ir 2 \E n 



/ &9.\f{r,6)\J l l 2 {r,e) J l ' 2 {r,e) 



drr 2 ( J dn\f(r,6)\J 1/2 (r,6) 



1/2 



snf (r) -4, rnr J Enl 



1/2 (223) 



where we use the Bishop-Gunther volume comparison theorem again, for J. Let us now 
make the observation that there are constants A+,A~, which depend only on 5 and k^s 
such that, 

(224) 



< — ^ < 



for r G [0,(5]. This can now be invoked at the second piece, giving us, 
5 



dr r e 



< 



/ dftl/M)!,/ 1 / 2 ^) J l l 2 (r,0) 



drr z / dtt\f(r,6)\'J(r, 



dn 



1/2 



^+(*:i,0) 



dre 



ft 2 c 2 



1/2 



<||/|| 2 (47r) 1 /^ + ( A;i ,0) 



2(m/n 2 C 2 ) 1/4 |^n| 1/4 



7.A 2 |fi n 



(225) 



If we choose A = (?i 2 i?/m) 1 ' / ' 3 |.E ri |~ 1 / 3 , the exponent goes to zero as \E n \ — > oo. For the 
first part of the integral we use the following characterization of essential supremum: let 
us define 

A(e) = n{{r G [0, A]| \r 3/2 F(r)\ > e}), (226) 



then we have 



Essup | r 3/2 F(r) | = inf{e|A(e) = 0} 

[o,A] 



(227) 



Let us use now F(r) = / g2 dfi/(r, 0), and using Bishop-Gunther bound for the first part 

as, 

-2, 







dr r 3 / 2 / dO |/(r,i 



n ' 2c 2 sn^ (r) 



,1/2 



(228) 
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which is smaller than; 



-2 

rM — 



^(fci.O) [ drr 3 / 2 [ dfi|/( 
Jo Js 2 

< A 2 + (k 1 ,0)( Essup\r 3/2 F(r)\) ( [ dr 

V [o,A] 1 \Jo 

< ^(fci,0)(Essup|r 3 / 2 F(r)| 

V [0,A] 



nr 2 | E„ 



-2 

e 



mr 2 |g n | 



r l/2 

1 



2(m/^C 2 ) 1 /4| J B n |i/4 



If we know take the limit A = {h 2 /mR 2 ) 1 / 3 ^^ 1 / 3 , we claim that the essential- 
suppremum goes to zero. To see this observe by Markov inequality [61] that 



A(e) < - [ A dr \r 3 ' 2 F( 
e Jo 



(\ ^l 2 /A / \ 2 \ ^/ 2 

J drr) I J drr 2 (J dO|/(r,0)| 

A 2 - \ X < 



;= / dr 



-sn 2 2 (r) / dO|/M)| 2 / dfi) 



£\/2\./o sn 2 2 (r) 1 
< —(tofPA+iO, k 2 ) I J* dr sn 2 2 (r) ^ dft |/(r, i 



< 7^( 4 -) 1/2 ^(°'^)^ A drr 2 ^ 

< ^(4^)1/2^(0,^)11/11,. 



dfl J(r,e)\f(r,< 



1/2 

v 1/2 
3M2 I 



hence for any e > 0, as A — > we can make A(e) = 0, thus the infimum goes to zero in 
this limit. As a result we see that the equation (1209|) is smaller than 

£^ {E ^ F(r)l) + ||/||2(4l) V 2 _^ e H &l ^/,,v. 



2{m/h?C 2 yl±\ [0 ,a] / 2(m//i 2 C 2 ) 1 /4 

2 — > as |i£ n | — > oo , 



e 

+ - 



S . , „. 1/2 (229) 



5 



9 



m d2( aj ,y)\E n \ 

dly e V 



and as a result goes to zero as desired. This completes the proof of the existence of the 
Hamiltonian in three dimensions. 
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